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Abstract. We calculate the ^'-meson energy spectrum in the T(1S) — > r/'ggg — > r/'X decay in the leading- 
order perturbative QCD in the static-quark limit for the orthoquarkonium. Our principal result is the 
extraction of parameters of the rj g* g effective vertex function (EVF) involving a virtual and a real gluon 
from the available data on the hard part of the ?/-meson energy spectrum. The perturbative-QCD based 
framework provides a good description of the available CLEO data, allowing one to constrain the lowest 

Gegenbauer coefficients B% and of the quark-antiquark and gluonic distribution amplitudes of the r\- 
meson. The resulting constraints are combined with the existing ones on these coefficients from an analysis 

of the rj — 7 transition form factor and the requirement of positivity of the EVF, yielding B^\^) = 

-0.008±0.054 and B { 2 9) 1 (jug) = 4.6±2.5 for til = 2 GeV 2 . This reduces significantly the current uncertainty 
on these coefficients. 



1 Introduction 

A quantitative description of the rare decays with the 
7/-meson production, such as B — > Tj'K^*\ B — » r)'X s , 
• T(1S) — ► vj X , requires an understanding of the rf g*g^*> 
\ effective vertex function (EVF), F n , g , g (*)(q 2 ,q 2 ,m 2 ,) [also 
called the rf — g transition form factor when one of the 
gluons is on the mass shell]. If one or both of the glu- 
, ons entering into the EVF in such decays are far from 
their mass shell, the dependence of the EVF on the gluon 
virtualities should be included in the theoretical anal- 
ysis of decays. For energetic ry'-meson, the QCD hard- 
scattering approach can be used to derive the required 
EVF 1,2,34]. As the fy'-meson has a relatively large mass, 
rrirj' — 958 MeV, its effect should be taken into account 
in applications of the t]'g*g^*' EVF to physical processes, 
in particular when the gluon virtualities are time-like g|. 
Moreover, the inclusion of the ry'-meson mass results in a 
pole-like form for the vj — g transition form factor: 

F n , g {p 2 ) = ^V rg (p 2 ,0X,) = _^L_±_Z (i) 

P m r l ' 

This form was introduced by Kagan and Petrov [o]. These 
authors also suggested to ignore the dependence of the 
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function H(p 2 ) on the gluon virtuality and approximate 
it by a constant value, H = 1.7 GeV -1 , resulting from 
the analysis of the J/ijj — > 7/7 decay |B] ■ Also it was subse- 
quently shown that the hard part of the r/-meson energy 
spectrum in the inclusive T(lS) — > rj'X decay jZ] is in 
a qualitative agreement with the spectrum measured re- 
cently by the CLEO collaboration Starting form the 
above observations, a quantitative analysis of the r/-meson 
energy spectrum in this decay was undertaken by us 9 
and its results are briefly summarized in this report. 



2 The //-Meson Wave-Function 

In the quark-mixing scheme, the Fock-state decomposition 
of the ry'-meson wave-function is as follows: 

W)=ew.<l>\rf q )+coB<f>\r/ a ) + \i/ g ), (2) 

where \rf q ) - \uu+dd)/V2, Ws) ~ \§s), and \rf g ) - \gg) are 
the light-quark, strange-quark and gluonic components, 
respectively. For an energetic ?/-meson in a process, its 
wave-function can be described in terms of the quark- 
antiquark, <j)f){x, Q 2 ), and gluonic, 4>^)(x, Q 2 ), light-cone 
distribution amplitudes (LCD As). In the LCD As above, 
x is the momentum fraction of one of the partons inside 
the meson and Q 2 is a typical hard scale of the process. 
Note that these LCDAs mix under the scale evolution. 

In most applications, approximate forms for the rf- 
meson LCDAs are usually employed in which only the 
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Pi = IV + *i - ?v 



-P/2,/3 



Fig. 1. A typical Feynman diagram describing the T(1S) 
999* (9* -> n'g) -> r/X decay 



first non-asymptotic term in both the quark-antiquark and 
gluonic component is kept: 

{x,Q 2 ) = 6xx [1 + 6(1 - 5xx) A 2 {Q 2 ) + . . .] , (3) 
$\x,Q 2 ) = 5x 2 x 2 (x-x)B 2 (Q 2 ) + ... . 

These LCDAs involve the Gegenbauer moments for which 
the following notation is used: 



A 2 (Q 2 



D. 



(</) 



a s (Q 2 



Jg) o(a) 

Pi B 2 



B 2 {Q 2 )= P fBf 



a s (Mo) 



a s (Q 2 ) 



B. 



(</) 



a s (Q 2 ) 



a s {Q 2 



(4) 
(5) 



The quantities 7+, 7^, p 2 q \ and p% are determined by 
the perturbative QCD while the Gegenbauer coefficients 
£?2 (^0) and Bif (Mo) are non-perturbative parameters. 
These coefficients have to be modeled or extracted from a 
phenomenological analysis of experimental data. 

The first attempt to estimate the Gegenbauer coeffi- 
cients B^lno) and B% (^o) was recently undertaken by 
Kroll and Passek-Kumericki 0. They performed a NLO 
theoretical analysis of the 7/ — 7 transition form factor and 
extracted the following values: 



B 2 } (t4 = 1 GcV ^) = 0.02 ±0.17, 
B 2 ] Oo = 1 Gev2 ) = 9 -° ±11-5, 



(6) 



from the CLEO [TUj and L3 [UJ data. This fit leaves an 
order of magnitude uncertainty on the coefficients. 

The inclusive T(1S) — > ?/X decay allows also to get 
an additional information on the Gegenbauer coefficients. 



3 Perturbative-QCD Analysis of T(1S) 
and Comparison with Data 



rfX 



One of the 18 diagrams describing the decay T(1S) — > 
999* (9* ~ * f/flO ~~ * m the leading order is presented 
in Fig. rjj the other 17 diagrams can be obtained from 
the above one by the permutations of the gluons in the 
intermediate (virtual) and final states. The static limit for 
the heavy quark and antiquark in the orthoquarkonium 
T(1S) state is used in the calculations. The total decay 
amplitude M.\T — > r/'ggg] is rather lengthy and can be 
found in j^j. 



Ho ~ 1.7 GeV- 
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Fig. 2. The ^'-meson energy spectrum in the T(1S) — > rj X 
decay 



rj' — 7 trans. 




T — ► rj'X decay 



-0.1 0.0 0.1 0.2 

Fig. 3. The resulting la contour (combined best fit), shown 
by the yellow (shaded) region, for the Gegenbauer coefficients 
estimated at the scale /ig = 2 GeV 2 from the data on the rj — 7 
transition form factor (solid curve) and T(1S) — > rj X decay 
(long-dashed and short-dashed curves) 



The yy'-meson energy spectrum can be theoretically de- 
termined as follows: 



1 1 



1 1 f <iki dk 2 dk.3 dpr/i 



(J) 



dn 

~cb ~ 3! (2tt) 8 2M J 2cji 2uo 2 2lo 3 2E v > 
x 6 {4) (V -h-k 2 -k 3 - p v >)6(z - 2E V ,/M) 



where is the three-gluon decay width of the Y(IS)- 
meson in the leading order: 



rl 0) = ^ (tt 



10 



9 



9)C F B F a 3 s (p 2 r ) 



M 2 



(8) 



Here, C F = {N 2 - 1)/(2JV C ), B F = (N 2 - 4)/(2JV c ), and 
/ir ~ M is a typical hard scale of the process. 
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Fig. 4. The ^'-meson quark-antiquark cj>jj, (x, fig) and gluonic cpjj/ (x, fig) LCD As as functions of x (left frame), and the resulting 
rj — g transition form factor (right frame). The solid curves correspond to the asymptotic quark-antiquark LCDA, while the 
LCDAs for the central values from the combined best-fit region of the Gegenbauer coefficients given in lilt are drawn as dotted 
curves. The LCDAs with the values B { 2 q) = 0.15 and B { 2 9) = 13.5, which are allowed within la from the analysis of the data 
on the rj — 7 transition form factor are presented as the dashed curves. Note that for this case the function F v i g (p 2 ) is no 
longer positive definite, as shown in the right frame 



As the low-z data are dominated by the fragmentation 
of gluons into the ?/-meson, following the CLEO analy- 
sis we concentrate on the last three (z > 0.7) and four 
(z > 0.6) experimental bins (see Fig.0). 

The fit of the B^\^q) and B^{j1q) parameters based 
on the last four experimental data points results in unac- 
ceptably large values of the minimum \ 2 [Hj. Thus, only 
the data in the last three bins with z > 0.7 are used in 
the analysis (quoted \ 2 corresponds to three degrees of 
freedom) and yield the following best fits — 2 GeV 2 ): 

B^Vu) = -0.89^, B^) = -2.86™ (9) 
i^Vo) = -1.09±i;^, ^ 9) (/^) = ll-53t 5 20 55 09 , (10) 

with x 2 — 2.45 and 2.37 for each set, respectively. 

The overlapping region in the [B% B% (a*o)] Pa- 

rameter space from the rj' — 7 transition form factor 3 
and the T(1S) — ► rj X decay 9 is presented in Fig. El The 
parameter space from the T(1S) — > rj'X decay fit is ob- 
tained by imposing the additional condition that F ri ' g {p 2 ) 
remains the positive-definite function for p 2 > to 2 , 9 
which is illustrated in Fig. 0] in term of both the LCDAs 
and the rj — g transition form factor. The resulting com- 
bined best fit of the Gegenbauer coefficients yields 0: 

^Vo) = -0.008±0.054, B^^l) =4.6 ±2.5. (11) 

Using instead the asymptotic LCDA yields \ 2 — 8.41; 
hence, the asymptotic LCDA is not favored by the current 
analysis. However, most of this \ 2 is contributed by a sin- 
gle experimental point (see Fig. As data on the hard 
part of the 7y'-meson energy spectrum are rather sparse, 
one can not exclude the asymptotic LCDA based on these 
data. Hopefully, experimental measurements will be im- 
proved soon to draw more quantitative conclusions. 



4 Summary 

The //-meson energy spectrum in the T(1S) — ► rfggg — ► 
r(X decay is calculated in the leading-order perturba- 
tive QCD in the static-quark limit for the T(15')-meson. 
The leading-twist (twist-two) quark-antiquark and gluonic 
LCDAs are used to describe the 7y'-meson wave-function. 
In the LCDAs, the asymptotic and the first non-asympto- 
tic terms are taken into account. An essential dependence 
of the energy spectrum on the Gegenbauer coefficients is 
observed. These Gegenbauer coefficients are determined in 
the large-z region (z > 0.7) of the f/-meson energy spec- 
trum from the recent CLEO data, however, the resulting 
la contour have a large dispersion. Combining this anal- 
ysis with the one of the rj — 7 transition form factor and 
requiring additionally that the EVF, F v i g (p 2 ), remains 
positive-definite in the entire p 2 > m 2 , region, yield much 
improved determination of the Gegenbauer coefficients. 

It remains to be seen if the so-determined rj g* g* EVF 
explains the data on the inclusive B — > rj X s decay. 

The work of A.Ya. P. has been supported by the Schweizerischer 
Nationalfonds. 
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